Zero- Turbulence Manifold in a Toroidal Plasma 
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Sheared toroidal flows can cause bifurcations to zero-turbulent-transport states in tokamak plas- 
mas. The maximum temperature gradients that can be reached are limited by subcritical turbulence 
driven by the parallel velocity gradient. Here it is shown that q/e (magnetic field pitch/inverse aspect 
ratio) is a critical control parameter for sheared tokamak turbulence. By reducing q/e, far higher 
temperature gradients can be achieved without triggering turbulence, in some instances comparable 
to those found experimentally in transport barriers. The zero-turbulence manifold is mapped out, 
in the zero- magnetic-shear limit, over the parameter space (7b, q/e, R/Lt), where 7_b is the per- 
pendicular flow shear and R/Lt is the normalised inverse temperature gradient scale. The extent 
to which it can be constructed from linear theory is discussed. 



Introduction. The heat loss that occurs as a result of 
turbulence driven by the ion temperature gradient (ITG) 
is one of the main obstacles to a successful fusion reactor. 
A large body of experimental work has demonstrated the 
effectiveness of strongly sheared equilibrium-scale flows 
in reducing this turbulence [l|-|3[. Numerical models [J- 
|7[ have demonstrated that by reducing the strength of 
the ITG instability, which drives the turbulence, and by 
shearing apart the turbulent structures, the radial gradi- 
ent of the flow component perpendicular to the magnetic 
field can indeed lead to a great reduction in the heat loss 
that results from a given temperature gradient. How- 
ever, Ref. y also demonstrated that the instability as- 
sociated with the parallel velocity gradient (PVG) [sHlOj 
could start to drive turbulence at higher flow gradients 
and prevent the complete suppression of the turbulent 
transport. 

More recent work has demonstrated that, at even 
higher flow shears, it is possible, at moderate temper- 
ature gradients, for the perpendicular velocity shear to 
overcome both the ITG and the PVG instabilities and 
completely quench the turbulence ll[. This quenching 



is most effective at zero magnetic shear |12Ml4J , a regime 
which has been associated in experiments with high con- 
finement of energy in the presence of strongly sheared 
flows [3|, ll5|. Refs. |12h14J also demonstrated the exis- 
tence, at zero magnetic shear, of a bifurcation to a high- 
temperature-gradient reduced-transport state, driven by 
a toroidal sheared flow. However, the maximum tempera- 
ture gradient that could be reached via such a bifurcation 
was found to be limited by the fact that turbulence was 
rekindled at high toroidal shear, in the form of subcriti- 
cal fluctuations driven by the PVG [9l4l3| . The question 
arises, which parameter regime is most favourable to the 
suppressing effect of the perpendicular flow shear and 
least favourable to the ITG and PVG drives? In other 



words, how can the temperature gradient which results 
from the transport bifurcation described in Refs. 12l4l4 
be maximised? 

At zero magnetic shear, the turbulence is subcritical 
for all nonzero values of the flow shear: there are no lin- 
early unstable eigcnmodcs, and sustained turbulence is 
the result of nonlinear interaction between linear modes 
which grow only transiently before decaying. A recent 
paper [10|, which studied this transient growth in slab 
geometry, demonstrated that at large velocity shears the 
maximal amplification exponent of a transiently grow- 
ing perturbation before it decays is proportional to the 
ratio of the PVG to the perpendicular flow shear. In a 
torus, this quantity is equal to the ratio of the toroidal to 
poloidal magnetic field components, or q/e, where q is the 
magnetic safety factor and e is the inverse aspect ratio. 
Therefore, if we conjecture that a certain minimum am- 
plification exponent is required for sustained turbulence, 
Ref. [10| predicts that there should be a value of q/e 
below which the PVG drive is rendered harmless. Below 
that value of q/e, it should be possible to maintain an 
arbitrarily high temperature gradient without triggering 
turbulent transport provided a high enough perpendicu- 
lar flow shear can be achieved. 

In this Letter, motivated by the possibility of reduced 
transport at low values of q/e, we use nonlinear gyroki- 
netic simulations to map out the zero-turbulence man- 
ifold, the surface in the parameter space that divides 
the regions where turbulent transport can and cannot 
be sustained. The parameter space wc consider is (7^, 
q/e, R/Lt), where 7^ is the normalised perpendicular 
flow shear: 7^ = u' /{q/e), where u' = dRto / dr / [vtui / R) 
is the toroidal shear, u the toroidal angular velocity, r 
the minor radius of the flux surface, Vtm the ion thermal 
speed and R the major radius, and where R/Lt is the in- 
verse temperature gradient scale length normalised to R. 



For brevity, we will refer to R/Lt as "the temperature 
gradient." We set the magnetic shear to zero, the regime 
we expect to be most amenable to turbulence quenching 
by shear flow [i[lll-[ll. 

We discover that reducing g/e is indeed uniformly ben- 
eflcial to maintaining high temperature gradients in a 
turbulence-free regime, and that values of R/ Lt can be 
achieved that are comparable to those ex per imentally ob- 
served for internal transport barriers [31, |l7| . 

In the next sections, having presented our numerical 
model and methodology, we will describe these results 
and discuss their physical underpinnings, as well as their 
implications for confinement in a toroidal plasma. Wc 
will show that linear theory of subcritical fluctuations 
[lOj can, with certain additional assumptions, provide 
good predictions of the nonlinear results. 

Numerical Model. To model the turbulence, we use 
the gyrokinetic equation [18| in the high-flow, low-Mach 
limit [l9| (i.e., the toroidal rotation velocity is ordered to 
be smaller than the sound speed but much larger than 
the diamagnetic velocity; Coriolis and centrifugal effects 
are neglected [35|, but velocity gradients are retained). 
We take the electrostatic limit and assume a modified 
Boltzmann electron response. The model used is iden- 
tical to that in Ref. [13|. The gyrokinetic system of 
equations is solved using the local nonlinear simulation 
code GS2 [Jol - E^ . As in Ref. [ll], we take the Cyclone 
Base Case parameter regime [23|, i.e., concentric circu- 
lar flux surfaces with e = 0.18, inverse ion density scale 
length R/Ln = 2.2 and ion to electron temperature ratio 
Ti/Te = l.|36| The magnetic shear is s = 0. The ratio q/e 
is varied by varying q alone. Collisions are included by 
means of a model collision operator, which includes scat- 
tering in both pitch angle and energy and which locally 
conserves energy, momentum and particles |24l . |25| . The 
resolution of all simulations was 128 x 128 x 40 x 28 x 8 
(poloidal, radial, parallel, pitch angle, energy). Note that 
relatively high parallel resolution was needed to resolve 
the PVG modes [H. 

Method. We wish to determine, in a three- 
dimensional parameter space (7^;, q/e, R/Lt), the 
boundary between the regions where turbulence can and 
cannot be sustained nonlinearly. We cover this space 
using four scans with constant q/e (Fig. [Ijb)), three 
scans with constant je (Fig- HIc)) and one scan with 
constant R/Lt (Fig. [ijd)). For each of these cases, wc 
consider multiple values of a second parameter and find 
the value of the third parameter corresponding to the 
zero-turbulence boundary. The boundary is defined as 
the point where both the turbulent heat flux and the 
turbulent momentum flux vanish. Thus, the location of 
each single point on the boundary is determined using 
on the order of ten nonlinear simulations. An example 
of this procedure is shown in Fig. [ija). In total, we per- 
formed more than 1500 simulations to produce the results 
reported below. 
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FIG. 1: (a) The simulations used to find the point on the 
manifold 7_e = 1.8, R/Lt = 15, q/e — 6.3, showing tlie heat 
flux vs. q/e at {je = 1.8, R/Lt ~ 15). The point on the 
manifold is the point where the heat flux drops to zero, (b-d) 
Sections through the critical manifold with parameters as in- 
dicated. Turbulence cannot be sustained for R/Lt < R/Ltc 
in (b,c), or for 7_b < ^ec in (d). The data points were found 
as illustrated in (a), and used to generate the manifold shown 
in Fig. H 



Because the turbulence that we are considering is sub- 
critical, there is always a danger that a simulation might 
fail to exhibit a turbulent stationary state because of an 
insufficient initial amplitude |26l . |27| . As we are not here 
concerned with the question of critical initial amplitudes 
we will consider a given set of parameters to correspond 
to a turbulent state if such a state can be sustained start- 
ing with a large enough perturbation. Therefore, all sim- 
ulations are initialised with high-amplitude noise. They 
are then run to saturation; close to the boundary a sim- 
ulation may need to run for up to t ~ 1000 R/vthi to 
achieve this. 

The critical curves obtained in this manner are plotted 
in Figs. Hfb-d). These curves, which effectively give the 
critical temperature gradient R/Ltc as a function of je 
and q/e, are then used to interpolate a surface, the zero- 
turbulence manifold, plotted in Fig. [5] The interpolation 
is carried out using radial basis functions with a linear 
kernel ^ (see also [H). 

Results. The results of the scan described above are 
displayed in Figs. [IJb-d). These three figures show, at 
fixed values of either 7^, R/Lt or q/e, the threshold in 
either R/Lt or q/e below which turbulence cannot be 
sustained; they arc. in effect, sections through the zero- 
turbulence manifold. 

Considering first Fig.[TJb), we sec that, at fixed q/e, the 
critical gradient R/Ltc first rises with je, as the perpen- 
dicular flow shear suppresses the ITG-driven turbulence. 




FIG. 2: The zero-turbulence manifold. Turbulence can be 
sustained at all points outside the manifold (that is, at all 
points with a higher temperature gradient and/or higher value 
of g/e than the nearest point on the manifold). This plot is 
made up from the sections shown in Fig. HJb-d) (heavy lines) 
and the manifold interpolated from them (thin grey mesh). 
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u' = dRu}/dr/{vt,^JR) 

FIG. 3: Contours of the zero-turbulence manifold plotted 
against the toroidal flow shear u = dRuj / dr / (vtM / R) = 
7_E/(5/e). The contours indicate the value R/Lt ~ R/Ltc 
below which turbulence is quenched. From top to bottom, 
the circles indicate approximate values of u' and q/e, corre- 
sponding to Ref. [li] (JET; R/Lt ~ 8), Ref. d (JET ITB; 
R/Lt ~ 17; note large discrepancy, see text) and Ref. [171 ] 
(MAST ITB; R/Lt - 10). 



and then falls — in most cases to — as the PVG starts 
to drive turbulence instead. This phenomenon was dis- 
cussed at length in |lll4l3| (indeed the curve for q/e = 8 
is taken from Il3|). Thus, for every q/e, there is an opti- 
mum value of the perpendicular flow shear ^e (and hence 
of the toroidal shear u') for which the critical tempera- 
ture gradient R/Ltc is maximised. We see that reducing 
q/e increases the maximum R/Ltc that can be achieved 
without igniting turbulence. Fig. [Ijc) shows that this 
rule applies for all considered values of flow shear j37| . 
This is to be expected, because lower q/e means weaker 
PVG relative to the perpendicular shear, allowing higher 
values of the perpendicular flow shear to suppress the 
ITG before the PVG drive takes over. 

Lastly, Fig.[IJd) shows the threshold in je above which 
the PVG can drive turbulence alone, without the help 



of the ITG; in other words, even configurations with a 
flat temperature profile would be unstable. At very high 
q/e, already a very small flow shear will drive turbulence; 
as q/e decreases, higher and higher values of je sue re- 
quired for the PVG turbulence to be sustained. It can- 
not be conclusively determined from this graph whether, 
as suggested by linear theory [lO|, there is a finite criti- 
cal value of q/e below which PVG turbulence cannot be 
sustained, i.e., a nonzero value of q/e corresponding to 
jEc -^ CO. However, for q/e ^ 7, the critical 7^ is far 
above what might be expected in an experiment |38| , and 
so the 7b — > cxD limit is somewhat academic. A definite 
conclusion we may draw is that at experimentally rele- 
vant values of shear, pure PVG-driven turbulence cannot 
be sustained for q/e ^ 7. 

The zero-turbulence manifold interpolated from the 
numerical data points is displayed in Fig. [2] The man- 
ifold comprises three main features: a "wall" where the 
critical temperature gradient increases dramatically at 
low q/e; a "spur" at low 7^, jutting out to high q/e 
(where, as 7^ increases, the ITG-driven turbulence is 
suppressed somewhat before the PVG drive becomes 
dominant), and finally the curve where the manifold in- 
tercepts the plane R/Lt = 0, whose shape is described 
above. 

Practical Implications and Comparison with Experi- 
ment. In order to illustrate better the implications of 
our findings for confinement, we plot, in Fig. [3j con- 
tours of R/ Ltc versus q/e and the toroidal flow shear 
u' — dRuj/dr/{vthi/R)- The basic message is clear: the 
lower the value of q/e, the higher the temperature gra- 
dient that can be achieved without igniting turbulence. 
Once we have obtained the lowest possible value of q/e, 
there is an optimum value of u' which will lead to that 
maximum R/Ltc- We note that the dependence of this 
optimum value of u' on q/e is not as strong as the de- 
pendence of the optimum value of 7^ on q/e (clearly this 
must be so because u' = {q/e)^E)- In a device with an 
optimised value of q/e, a near maximum critical temper- 
ature gradient would be achievable for u' ^5, shears 
comparable to those observed in experiment [sl, ll6|, ll7[ . 

While simulation results obtained for Cyclone Base 
Case parameters are not suitable for detailed quantita- 
tive comparison with real tokamaks, it is appropriate to 
ask whether our results are at all compatible with experi- 
mental evidence. For an internal transport barrier (ITB) 
in MAST, Ref. [I^ reports R/Lt ~ 10 at q/e - 4.6 ^ 
and u' ^ 2.4. This is comparable to the critical values 
shown in Fig. H In JET, Ref. [H reports R/Lt ~ 8 at 
q/e ^ 11 and u' ^ 3.6, again reasonably close to what we 
would have predicted. However, an ITB in JET studied 
by Ref. [3 achieved R/Lt -- 17 at q/e - 10 and u' - 4.1 
— substantially higher than our R/Ltc at the same val- 
ues of q/e and u' . Note, however, that Ref. [3[ reports 
that the shear was dominated by an enhanced poloidal 
flow, an effect which is not included in our numerical 



model. 

Relation to Linear Theory. Since the mapping of the 
zero-turbulence manifold using nonlinear simulations is 
computationally expensive, wc may ask whether linear 
theory can predict marginal stability. The question is 
also interesting in terms of our theoretical understanding 
of subcritical plasma turbulence. It is clear that in a situ- 
ation where perturbations grow only transiently, existing 
methods based on looking for marginal stability of the 
fastest growing eigenmode will not be applicable. In Ref. 
[lOJ I , we considered these transiently growing modes in a 
sheared slab, and posited a new measure of the vigour of 
the transient growth: A'max, the maximal amplification 
exponent, defined as the number of e- foldings of transient 
growth a perturbation experiences before starting to de- 
cay, maximised over all wavenumbers. It appears intu- 
itively clear that in order for turbulence to be sustained, 
transient perturbations must interact nonlinearly before 
they start to decay. We may then assume that a satu- 
rated turbulent state will exist if iV^ax ^ Nc, where Nc is 
some threshold value of order unity. The zero-turbulence 
manifold is then the surface Ny^axijE, g/e, R/Lt) = Nc- 

We now test this idea by calculating A'^max for lin- 
ear ITG-PVG-driven transient perturbations in a slab, 
using the code AstroGK [29| to solve the linearised gy- 
rokinetic equation, as done in Ref. [Ifll. Fig. ^&) 
shows that for each value of 7^; and a range of q/e, 
it is possible to choose Nc{'Je) such that the equation 
.^max(7B, g/e, R/Lt) = N(. correctly reproduces the crit- 
ical curve R/Ltc{iI^) obtained as a section of the zero- 
turbulence manifold at that value of 7^- However, Nc 
does have a strong dependence on ^e, shown in Fig.|4l[b), 
ranging from Nc ^ 0.5 at 7^ > 2 to iVc ^ oo as 7^ — >■ 
(the latter is an expected result: at 7^ = 0, there is a 
growing eigenmode, so either iVmax = cx) or there is no 
growth at all) . It is not clear if Nc tends to a finite limit 
as 7^; ^- 00, but, similarly to the existence of a critical 
value of q/e as ^e -^ 00, this is a somewhat academic 
question because such a limit would be achieved (or not) 
at ^E too large to be experimentally achievable. 

The practical conclusion of this exercise is that all that 
appears to be required to determine the two-dimensional 
dependence of R/Ltc on 7^ and q/e is finding Nd^fE) 
using a nonlinear scan at a single value of q/e; thus, the 
number of parameters in the nonlinear scan is reduced 
by one. 

Discussion. We have presented two key results. 
Firstly, and principally, we have calculated the shape of 
the zcro-turbulcncc manifold, the surface that divides the 
regions in the parameter space (7^, q/e, R/Lt) where 
subcritical turbulence can and cannot be nonlinearly sus- 
tained. We have described the shape of this manifold 
and its physical origins, and presented its two implica- 
tions for confinement in toroidal plasmas: that reducing 
the ratio q/e, i.e., increasing the ratio of the poloidal to 
the toroidal magnetic field, improves confinement at ev- 
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FIG. 4: (a) The critical temperature gradient R/Ltc vs. q/e 
for different values of 7_e, showing both R/Ltc obtained from 
the interpolated manifold, and R/Ltc such that A'^niax = Nc, 
with Nc suitably chosen for each 7_b, as shown in (b). 



ery nonzero value of je, and that at fixed q/e, there is 
an optimum value of 7e (that is, an optimum value of 
the toroidal flow shear u' — dRoj/dr/{vthi/R)) at which 
the critical temperature gradient is maximised, in some 
instances to values comparable to those observed in in- 
ternal transport barriers [3|, ll7|. How to calculate the 
heat and momentum fluxes that would need to be in- 
jected in order for such optimal temperature gradients 



to be achieved was discussed in Ref. 14 



Secondly, we have shown that the zero- 
turbulence manifold can be parameterised as 
Nma^'-fE, q/e, R/Lt) = Nc{je), where A^max is the 
maximal amplification exponent of linear transient per- 
turbations (calculated from linear theory) and Nc must 
be fit to the data. Thus, using a single scan at constant 
q/e to determine Nc^je) appears to be sufficient for 
calculating the full two-parameter dependence of the 
critical temperature gradient. Obviously, the need to fit 
NdjE) indicates a limitation of our current theoretical 
understanding of the criterion for sustaining subcritical 
turbulence in a sheared toroidal plasma. The results 
reported here provide an empirical constraint on future 
theoretical investigations. 

Another avenue for future investigations is determin- 
ing the dependence of the zero turbulence boundary on 
some of the parameters that were held fixed in this work: 
Ti/Tc, magnetic shear, and, more generally, the shape of 
the flux surfaces, density gradient, inverse aspect ratio e 
(separately from q), etc. Mapping out the dependence 
just on ^E, g/e and R/Lt took approximately 1500 non- 
linear simulations at a total cost of around 4.5 million 
core hours. Adding even two or three more parameters 
to the search would take computing requirements beyond 
the limit of resources today, but not of the near future. 
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A temperature ratio of 1 is appropriate for both lower 
power and future reactor-like conditions, but not for 
high-performance shots in current devices [la . |33 . |33| . 
The increase at 7_b = cannot, of course, be due to re- 
duction of the PVG; we assume that this occurs because 
of the simultaneous reduction of the maximum parallel 
length scale qR in the system, leading to weaker ITG 
turbulence; see [34| . 

By order of magnitude, je ~ M/q, where M is the Mach 
number of the toroidal flow. Thus, values of 7_e much 
above unity are unlikely to be possible. 
The ratio of toroidal to poloidal field in MAST can be 
smaller on the outboard side, so the effective value of 
q/e for locating this case on the zero-turbulence manifold 
might be smaller than quoted. 



